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ABSTRACT 

We use the extensive catalog of dark matter haloes from the Millennium simulation to 
investigate the statistics of the mass accretion histories (MAHs) and accretion rates 
of ~ 500,000 haloes from redshift z = to 6. We find only about 25% of the haloes 
to have MAHs that are well described by a 1-parameter exponential form. For the 
rest of the haloes, between 20% (Milky- Way mass) to 50% (cluster mass) experience 
late-time growth that is steeper than an exponential, whereas the remaining haloes 
show plateau-ed late-time growth that is shallower than an exponential. The haloes 
with slower late-time growth tend to reside in denser environments, suggesting that 
either tidal stripping or the "hotter" dynamics are suppressing the accretion rate of 
dark matter onto these haloes. These deviations from exponential growth are well 
fit by introducing a second parameter: M{z) oc (1 -I- z)^e~^^. The full distribution 
of P and 7 as a function of halo mass is provided. From the analytic form of M(z), 
we obtain a simple formula for the mean accretion rate of dark matter, M, as a 
function of redshift and mass. At z = 0, this rate is 42M0yr~^ for IO^'^^Mq haloes, 
which corresponds to a mean baryon accretion rate of = 7M0yr~^. This mean 



rate increases approximately as (1 



\1.5 



at low z and (1 -I- z) ° at high z, reaching 



Mb = 27,69, and 140 M0yr ^ at z = 1, 2, and 3. The specific rate depends on halo 
mass weakly: M/M oc M^ -^^"^. Results for the broad distributions about the mean 
rates are also discussed. 



1 INTRODUCTION 



The mass growth history is a basic property of dark matter 
haloes. Haloes in numerical simulations are seen to be assem- 
bled through a number of processes: mergers with compa- 
rable mass haloes ("major mergers"), mergers with smaller 
satellite haloes ("minor mergers"), and accretion of non-halo 
material that is composed of either haloes below the numeri- 
cal resolution or diffuse particles. Following the mass history 
of the most massive progenitor halo as a function of redshift 
z is a useful way to quantify a halo's mass assembly history. 
These mass accretion histories (MAHs, or M{z)) are impor- 
tant for statistical studies of the distributions of halo forma- 
tion redshifts, and the correlations between formation time 
and other halo properties such as environment, concentra- 
tion, substructure fraction, spin, and relative contributions 
to mass growth from major vs minor mergers. Moreover, the 
time derivative of the MAH gives the mass growth rate of 
dark matter haloes, which is directly related to the accre- 
tion rate of baryons from the cosmic web onto dark matter 
haloes. 

A number of earlier papers have investigated various as- 



were seen, with late-forming haloes being less concentrated. 
The scatter in c was attributed to the scatter in z/. An al- 
ternative 2-parameter fitting function was demonstrated by 
van den Bosch ( 2002 1 to be superior to a 1-parameter fit to 



haloes in a simulation with the same particle number in a 
141h~^ Mpc box. 

The relationship between halo structure and accretion 



was further addressed in Zhao et al. ( 2003 \ and Zhao et al. 
(20031, where the redshift dependence of c was observed 



to be more complicated than a simple proportionality. |Ta- 



sitsiomi et al. (20041 examined 14 haloes. 



ranging m mass 
^*/i-^Mr) and 



from group to cluster scale (.58 to 2.5 x 10 n iviq 
also found that a 2-parameter fit for Af(z) worked better. 
Cohnfc White] (2005) studied the mass accretion histories of 



1500 cluster-sized haloes and characterized several prop- 
erties of galaxy cluster formation. 



Maulbetsch et al. ( 2007 1 studied the environmental de- 
pendence of the formation of ~ 4700 galaxy-sized haloes 
(above lO"/i"^Af0) in a 50 /i ^ Mpc simulation box. In 
higher-density environments, they found the haloes to form 
earlier with a higher fraction of their final mass gained via 
major mergers. |Li et al." 12008) studied 8 different definitions 



pects of the halo MAHs. For instance, ^Wechsler et al. ( 2002 1 
analyzed ~ 900 haloes (above 10^^h~ Mq at z = 0) in a 
ACDM simulation in a 60 h'^ Mpc box with 256^ parti- 
cles. The values from a 1-parameter fitting function for the 
MAHs were presented for 8 haloes. Clear correlations be- 
tween the formation redshift z/ and concentration c of haloes 



of halo formation time using the haloes from the Millen- 
nium simulation ( [Springel et al.|2005[ ). The motivation was 
to search for halo formation definitions that better charac- 
terize the downsizing trend in star formation histories, as 
opposed to the hierarchical growth of haloes in the ACDM 



cosmology. Zhao et al. (2008j) (Z08) investigated the mean 
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MAH in different cosmological models - scale-free, ACDM, 
standard CDM, and open CDM - and searched for scaled 
mass and redshift variables that would lead to a universal 
fitting form for the median MAH for all models. 

The results in these earlier papers were presented either 
for M{z) of a handful of individual haloes, or for the global 
mean growth of a selection of haloes. Our aim here is to 
quantify systematically the diversity of growth histories and 
rates using the ~ 500, 000 z = haloes with M > W^^Mq 
(i.e. above 1000 simulation particles) and their progenitors in 
the Millennium simulation. Over this large range of haloes, 
we find that an exponential fit does not adequately capture 
the behavior of halo growth. Many haloes experience large 
changes in the rate at which they accrete mass. Some haloes 
grow more slowly at late times, and occasionally even lose 
mass, while other haloes undergo late bursts of growth. All 
of these MAHs are poorly fit by an exponential, and suggest 
the need for a fitting form with more fiexibility. We find it 
helpful to classify the MAHs into four types based on their 
late-time accretion rate. The large ensemble of haloes allows 
us to quantify the mean values as well as the dispersions of 
the mass accretion rates and halo formation redshifts as a 
function of mass and redshift. 

This paper is organized as follows. Sec. 2 provides some 
background information about the haloes in the Millen- 
nium simulation and describes how we construct halo merger 
trees. This post-processing of the Millennium public data is 
necessary for identifying the thickest branch (i.e. the most 
massive progenitor) along each final halo's past history. The 
masses of these progenitors will then allow us to quantify the 
MAH, M{z). In Sec. 3, we first assess the accuracy of the 
1-parameter exponential form for AI{z). We then propose a 
more accurate two-parameter function for M{z) and classify 
the diverse assembly histories into four broad types accord- 
ing to their late-time growth behavior. We further quantify 
the statistics of the two fitting parameters, providing (in the 
Appendix) algebraic fits for their joint distributions that can 
be used to generate Monte Carlo realizations of an ensemble 
of halo growth tracks. The applicability of M{z), which is 
derived for z — haloes, for the mass accretion history of 
higher-redshift haloes is discussed in Sec. 3.3. Sec. 4 is fo- 
cused on the statistics of the mass accretion rates. A simple 
analytic expression is obtained for the mean accretion rate, 
^M^, of dark matter as a function of halo mass and red- 
shift. The dispersions about the mean rates are significant, 
as evidenced by the differential and cumulative distributions 
of M presented here. Sec. 5 discusses the mean and the dis- 
tribution of the halo formation redshift as a function of halo 
mass. In Sec. 6 we report the correlations of MAHs with 
halo environment, the last major merger redshift, and the 
fraction of haloes' final masses assembled via different types 
of mergers. 



HALO MERGER TREES IN THE 
MILLENNIUM SIMULATION 



The Millennium simulation ( [Springel et al.|2005[ ) provides a 
database for the evolution of roughly 2x10' z = dark mat- 
ter haloes from redshifts as high as z = 127 in a 500h~^ Mpc 
box usmg 2160^ particles of mass 1.2 x IO^'Mq (all masses 
quoted in this paper are in units of Mq and not H'^Mq). 



It assumes a ACDM model with fi™ = 0.25, Qb = 0.045, 
Qa = 0.75, h = 0.73, and a spectral index of n = 1 for the 
density perturbation power spectrum with a normalization 
of (78 = 0.9. 

Dark matter haloes are identified with a friends-of- 
friends (FOF) group finder ( Davis et al.||198"5 l with a link- 
ing length of 6 = 0.2. Throughout this paper we use the 
number of particles linked by the FOF finder to define the 
halo's mass. Once identified, each FOF halo is then bro- 
ken into gravitationally bound substructures (subhaloes) by 
the SUBFIND algorithm (see [Springel et al.|[200T| ). These 
subhaloes are connected across the 64 available redshift out- 
puts: a subhalo is the descendant of a subhalo at the preced- 
ing output if it hosts the largest number of the progenitor's 
bound particles. The resulting subhalo merger tree can be 
used to construct merger trees of FOF haloes, although we 



have discussed at length in Fakhouri & Ma ( 2008 2009 1 the 



complications due to halo fragmentation and have presented 
comparisons of several post-processing algorithms that han- 
dle fragmentation events. 

Our results in this paper are based on the stitch-3 post- 



processing algorithm described in Fakhouri & Ma ( 2008 1 



In this algorithm, fragmented haloes that remerge within 3 
outputs after fragmentation are stitched into a single FOF 
descendant; those that do not remerge within 3 outputs 
are snipped and become orphan haloes. After applying the 
stitch-3 algorithm, we extract the mass accretion history, 
M{z), of each halo at z = (or at any higher redshift) by 
following the halo's main branch of progenitors. We have 
compared the resulting AI{z) and formation redshifts to 
those obtained from the alternative algorithms (e.g., "snip," 
"split," and subhalo vs FOF mass) discussed in Fakhouri & 



Ma ( 2008 2009 1. We find the systematic variations to all be 



within 5-10% of the stitch-3 values of these quantities. 



3 FITTING MASS ACCRETION HISTORIES 
3.1 Previous MAH Forms 

To quantify the limitations of the exponential fit in cap- 
turing halo growth, consider the formation redshift Zf, here 
defined as the redshift at which M{z) is equal to Mo/2. 

For the 1-parameter exponential form (e.g. |Wechsler| 
et al.|2002| ) 



M{z) = Moe""", 
the parameter a is simply related to by 



ln(2) 



(1) 



(2) 



We have compared 2/ as determined by the exponential fit to 
each halo's M{z) from the simulation with the Zf determined 
directly from the M{z) tracks such that M{zf ) = M(0)/2 
(using interpolation between output redshifts). We find the 
exponential fit to err systematically in its determination of 
Zf, significantly overestimating the formation redshift for 
haloes that form recently and underestimating it for haloes 
that form early. The mean value of 2/ from the exponential 
fit, for instance, is 0.3 higher than the actual value for young 
haloes and is 0.8 lower for old haloes across all masses. 

A more complicated functional representation of MAHs 
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was put forth by van den Bosch ( 2002 I 



log 



M(z) 
Mo 



= -0.301 



log(l + 
log(l + ^/) 



(3) 



where u and Zf are fitting parameters. The use of an ad- 
ditional parameter provided significant improvement in the 
quality of the fits for many MAHs, especially those that 
formed late. This 2-parameter form, however, is not flexi- 
ble enough to handle haloes that have lost mass, as it can- 
not take on values that would give M{z)/Mq greater than 
1. Moreover, over the sample of haloes tested in [van deri] 
Bosch ( 2002 I, comparison between the goodness-of-fit of this 



two-parameter form and the exponential fit showed that the 
exponential fit actually performs better for early forming 
haloes. 




0.06 0.08 0.10 
RMS Deviation 



3.2 A Revised MAH Form 

To address the need for a fit that is both effective and simple, 
we find a 2-parameter function of the form 



M{z) = Mo{l + zfe~ 



(4) 



to be versatile enough to accurately capture the main fea- 
tures of most MAIfs in the Millennium Simulation. This 



form has also been studied in Tasitsiomi et al. (20041 for 



cluster-mass haloes, but it has not been tested over a large 
number of haloes of different mass. The form reduces to an 
exponential when /3 is 0, and 7 in this case is simply the 
inverse of the formation redshift: 7 = ln{2)/zf. A large frac- 
tion of the haloes, however, are better fit when the additional 
factor of (1 -I- z)'^ in equation Q is included. In general, /3 
can be either positive or negative, but 7 0. We find the 
combination /3 — 7 to be a useful parameter for character- 
izing these MAHs as /3 — 7 gives the mass growth rate at 
small redshifts: 



dln(M(z)) _ 13 



dz 



- 7 /3 - 7 + C(z) . 



(5) 



This late-time trend can be used to characterize the MAH 
as described below. 

To obtain the best-fit values for /3 and 7 in equation Q, 
we have performed a x^-like minimization of the quantity 



N ^ 



[M{z,)/Mo-{l + Zif 
M{z,)/Mo 



(6) 



where the sum is over the A'^-available simulation redshift 
outputs at Zi{i = 1, N) for each halo. The choice of the 
factor M{zi) in the denominator is akin to assuming Pois- 
sonian errors for halo masses. We found this choice to be a 
suitable middle ground between minimizing simply the sum 
of squares and minimizing the fractional deviation (i.e. with 
a factor of M'^{z) in the denominator). The former tended to 
fit the finely sampled low- 2: points well at the expense of the 
sparsely spaced high-z points, whereas the latter tended to 
do the opposite. Equation (|6|, on other hand, provides rea- 
sonable fits for the entire history of the halo growth. Fig. [l] 
shows the cumulative distribution for the rms deviation of 
the fits from the Nbody data (normalized by Mq) for all 
~ 500, 000 2 = haloes. The deviation is less than 6% for 
over 75% of the haloes, and only a few percent of haloes have 
deviations larger than 10%. Of this most poorly fit subset 
of haloes, nearly half underwent mass loss at late times. 



Figure 1. Cumulative distribution for the RMS deviations be- 
tween the fitting formula eq. |4]l and the Millennium Nbody out- 
put for the mass growths of 478,781 2 = haloes (solid black) and 
the four sub- types (colored curves) listed in Table 1. The figure 
shows that the fits perform well overall: about 75% of the haloes 
have RMS deviations less than ~ 6%. 



Type 



Criteria 



Characteristics 



Xl 



I 

II 
III 

IV 



< 0.35 
/3 - 7 < -0.45 
-0.45 < /3 - 7 < 
/3-7 > 



Good exponential 1.09 

Steep late growth 1.61 

Shallow late growth 2.15 

Late plateau/decline 3.31 



Table 1. MAHs are categorized based upon the best fit param- 
eters /3 and 7 of equation Q . Categorization is done in order by 
type; thus MAHs that satisfy the criteria for both Type I and 
Type II belong to Type I. The right-most column is the mean 
of Xi/xi' ratio of the computed for the 1-parameter ex- 
ponential form to the x^ computed for the 2-parameter form in 
equation |4jl. Values > 1 imply that the 2-parameter form pro- 
vides a more accurate fit than the 1-parameter form. 



As expected, the fits become progressively worse at higher 
redshifts; for over 75% of haloes, the maximum fractional 
deviation between the fits and Nbody results occurs above 

2 = 4. 

We suggest that the parameters (5 and 7 allow for rough 
classifications of MAHs into a few basic groups, summarized 
in Table [l] The classification scheme is quite straightfor- 
ward. Fits with small values for /3 indicate a weak contribu- 
tion from the (1 4- 2) term, and deviate minimally from an 
exponential curve. These haloes with < 0.35 are labeled 
Type I. 

The rest of the classifications are dependent upon the 
value of /3 — 7. The motivation for this is the fact that the 
difference represents the value of the derivative at 2 = 0, 
as noted in equation (|5|. Hence Type II haloes, defined to 
be those haloes with /3 — 7 < —0.45, feature steep growth 
at late times, typically steeper than can be captured by an 
exponential fit. 

Type III haloes have fit parameters that fall in the range 
—0.45 < /3 — 7 < and exhibit flat late time growth. Like 
Type II, these tend to deviate from the fit that would be 
found using the exponential form, but Type HI haloes do so 
in the opposite direction to Type II haloes. A typical Type 
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Figure 2. Average mass accretion histories (top) and their 
derivatives (bottom) for the total population (black solid) and 
the four types of halo growths listed in Table ^ The average is 
taken over the M{z) for galaxy-sized haloes with masses between 
2.1 X IO^^Mq and 3.3 X IO^^Mq in the Millennium simulation, al- 
though haloes at different mass show similar behavior. In the top 
panel, the set of curves with the lighter shading shows the average 
M[z) computed from our fits of equation Q to each halo's MAH. 
The bottom panel illustrates that the late-type growth rates dif- 
fer greatly, ranging from d\r\M/dz ~ 1.2 for Type II to 0.1 for 
Type IV. 



Ill halo has undergone limited growth during recent times, 
sometimes after a spurt of growth at earlier times. 

Type IV, with /3 — 7 > 0, represents the most extreme 
deviation from an exponential. The majority of Type IV 
haloes have shed mass, some of them by significant amounts. 
Some Type IV haloes have merely seen their growth slow 
down like the Type III haloes, but over a more significant 
period of time. As such. Type IV haloes are extreme cases 
of Type III haloes, perhaps representing the future growth 
for some Type III haloes. 

The boundaries delineating these classifications are 
rough guidelines at best. For example, consider the defini- 
tion for Type I of |/3| < 0.35. For the largest values of /3 
in this group, which should be considered the worst of the 
"good exponentials" that constitute Type I, the fractional 
difference between the formation redshift as determined by 
the simple exponential and the modified exponential is a lit- 
tle under 8%. The agreement is not perfect, but the two fits 
are similar enough for these haloes that the use of the power 



Mass Range 


Halo Number 


Type I 


TT 
ii 


TTT 
ill 


IV 


1.2 to 2.1 


191421 


29% 


27% 


32% 


12% 


2.1 to 4.5 


143356 


27% 


29% 


32% 


12% 


4.5 to 14 


95744 


24% 


34% 


31% 


11% 


14 to 110 


43089 


20% 


42% 


26% 


11% 


> 110 


4787 


18% 


57% 


17% 


8% 




478781 


27% 


31% 


31% 


11% 



Table 2. Within each mass range, the percentage of haloes that 
belong to each type are provided. For each type, there is a notice- 
able trend with mass, though the strength of the trend varies. 



law parameter adds little. Of course, there is no reason why 
we should not instead demand that the formation redshifts 
differ on average by no more than 5%, or perhaps 10%. In 
the end, the combination of the formation redshift metric 
and a couple of others for comparing the fits suggested that 
demanding < 0.35 was inclusive enough to capture the 
majority of haloes for which an exponential is an adequate 
fit, without unduly diminishing the integrity of the group. 

Fig. [2] compares the shapes of the average MAH for 
haloes of galaxy-size mass from the Millennium simulation 
for the overall distribution and for each type. The bottom 
panel shows the derivative d In M/dz to highlight the differ- 
ent late-time accretion rates among the four types. Haloes 
of other mass show similar behavior. Clearly, the late time 
growth rate is an important factor in distinguishing haloes 
from one another. The average MAH for Type I haloes is 
quite similar to that of the overall distribution, which indi- 
cates that the average MAH is approximately exponential. 
However, the behavior of about 75% of individual haloes 
deviates from an exponential noticeably. This fact is quan- 
tified in the right- most column of Table [l] where the ratio 
of for the exponential fit to the 2-parameter fit is seen to 
increase with the MAH types. 

Since the mean MAH is approximately exponential, the 
accretion rate din M/dz averaged over the whole population 
is also nearly independent of redshift (black solid curves 
in Fig. [2| when expressed in units of per redshift, with 
d\nM/dz being between 0.6 and 0.7 for 2 = up to 5. This 
weak dependence on redshift is similar to that of the halo 
merger rates (per unit z) reported in |Fakhou ri fc Ma'(2008[). 
The different types of haloes, however, show significant dis- 
persions in the late-time accretion rates, with d In M/ dz be- 
ing as high as 1.2 for Type II and as low as 0.1 for Type IV 
at 2 « 0. 

For each of the mean profiles shown in Fig. [2] we have 
fit the analytic form in equation Q. The best-fit values of 
(/?, 7) are (0.10, 0.69) for all haloes, and (-0.04, 0.54), (-0.9, 
0.35), (0.62, 0.88), and (1.42, 1.39) for each of the four types, 
respectively. 

The statistics of the 478,781 2 = haloes (above 1000 
particles, or a mass of 1.2 X lO^^M©) belonging to each MAH 
type across different mass bins are given in Table |2] and 
Fig. [3] Exponential MAH (Type I) is seen to apply to only 
20 to 30% of the haloes. There is also interesting dependence 
of the type on halo mass. Most notably. Type II haloes fea- 
ture a strong dependence on mass, where the fraction rises 
from 27% at ~ 10^^ M© to 60% at > IO^^Mq. Cluster-size 
haloes therefore not only form late, which is a natural con- 
sequence of the ACDM cosmology, but the majority of their 
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Mo (Ms) 



Figure 3. Cumulative fraction of lialoes belonging to a given 
MAH type as a function of halo mass (e.g. the magenta type-IV 
curve includes the contributions from the three other types) . The 
exponential form (Type I) is a good fit for only 20% to 30% of 
the haloes at all masses. The type II fraction shows a strong mass 
dependence, reaching ^ 60% for cluster-mass haloes. 



mass accretion rates is also faster than an exponential at 
low redshifts. 

The results presented thus far are for the mean MAH 
and mean values of (/9, 7). We find, however, significant dis- 
persions about the mean behavior that are also important 
to characterize. For completeness, we show the distributions 
of our best-fit (/?, 7) for all halo MAHs in the Appendix and 
Fig. |A1[ We also present there an accurate fitting form that 
we have obtained for the two-dimensional probability distri- 
bution of (3 and 7 as a function of halo mass. This formula 
can be used to generate a Monte Carlo ensemble of realistic 
halo growth histories. The details of the formula, its usage, 
and comparison to the Millennium data are described in the 
Appendix. We emphasize that the results presented for the 
rest of this paper are obtained from the Millennium haloes 
directly rather than from this Monte Carlo realization. 

3.3 MAHs for Haloes at Higher Redshifts 

The MAHs presented thus far are obtained from the main 
branches of the descendant haloes at z = 0. Thus, for a 
higher redshift zi > 0, the distribution of M(2i) contains 
only information about the main branch progenitors, which 
is a subset of all the haloes at zi since many haloes do not 
belong to main branches. 

Since the formation of higher-redshift galaxies and their 
host haloes is of much interest, it is useful to quantify the 
behavior of MAHs for haloes at zi , where zi > 0. In partic- 
ular, we ask whether the mean MAH for haloes of mass Mi 
at zi for z > zi can be related to the MAHs of haloes at 
z = that we have studied thus far. 

We find that the mean MAH of haloes of mass Mi at zi 
is nearly identical to the mean MAH of haloes at z = that 
satisfy M{zi) = Mi. That is, the mean MAH for z > zi oi 
the main branch subset with mean mass M{zi) = Mi at z\ is 
very similar to the mean MAH of the complete population of 
haloes with mean mass Mi at 2:1 . As a specific example, the 
mean MAH of the 10^^ M© z — Q haloes in the simulations 
had the value M{z = 1) = 4.5 x lO^^M© at 2 = 1. We find 




Figure 4. For haloes of mass Mi at redshift zi , the y-axis plots 
their corresponding mean mass Mq today. Four values of z\ are 
shown: (solid), 0.5 (dashed), 1 (dot dashed), and 2 (dotted). 
This mapping allows one to use equation (j4]| for the MAH of 
higher-redshift haloes (see text). 



that the mean MAH of these lO^^M© haloes at z > 1 is 
nearly identical (within 2%) to the 2 > 1 evolution of the 
mean MAH of all Mi = 4.5 x lO^^M© haloes at 2 = 1. This 
property for the mean MAH is in fact a natural consequence 
of the Markovian nature of the Extended Press-Schechter 
theory (see, e.g.. Sec 2.3 of |White||1994 1. 



This self-similar property implies that in order to study 
the MAH properties of haloes with mass Mi at redshift 21 , 
one simply needs to determine which set of haloes at 2 = 
have M(2i) = Mi. In particular, one needs to compute the 
average mass Mo of the haloes at 2 = that map onto 
M(2i) = Ml at 2i. This mapping is shown in Fig. |4] with Ml 
along the x-axis and Mo along the y-axis for 21 = 0, 0.5, 1, 
and 2. Note that Mo = Mi at 21 = by construction, and 
as 21 increases, the mass Mo that maps onto some fixed Mi 
by redshift 21 also increases. 

We note that the mapping in Fig. [4] implies that haloes 
of some mass Mi at some redshift 21 > do not have the 
same shape of MAH as haloes of mass A/o = Mi at 2 = 0. 
That is, the MAH of a IO^Mq halo at 2 = and the MAH 
of a 10^"^ M0 halo at 21 > are not simply related by a 
shift from 2 to 2 — 21 in equation Q. This is because haloes 
at higher 21 have a relative formation redshift 2/ — 21 that 
is smaller than haloes of the same mass at 2 = 0. This 
result is not surprising since haloes of the same mass at 
different redshifts in the ACDM model represent different 
part of the mass spectrum and are not generally expected 
to have identical properties. 

We have tested the self-similar property of the fitting 
form of Z08 (using their online code) by comparing their 
mean MAH for 2 = lO^M© haloes and the MAH for their 
M(2 = 1) haloes at 2 > 1. Their latter MAH is higher than 
the former by about 15%, while ours differ by less than 2%. 



4 MASS ACCRETION RATES: MEAN AND 
DISPERSION 

Having quantified M[z) in Sec. 3, we now examine its time 
derivative - the mass accretion rate - in more detail. In 
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Figure 6. Differential (top) and eumulative (bottom) distributions of the accretion rates of cosmic baryons, A/;,, for four halo masses 
(left to right). Within each panel, the accretion rates at four redshifts z = (solid), 0.5 (dashed), 1 (dashed dotted), and 2 (dotted) 
are shown, where the distributions are seen to broaden significantly with increasing z. The vertical axis labels the number of haloes per 
comoving Mpc^ at or above a given Af;,. 



particular, we would like to obtain a general formula for 
the mean accretion rates of dark matter for a wide range of 
halo mass and redshift. To achieve this, we note that our 
analytical form in equation Q for individual halo MAHs 
gives: 



M 
M 



= 0.10/iGyr"' [7(1 + z)-(]] y^Qmil + + (7) 



where Qm and f^A are the present-day density parameters 
in matter and the cosmological constant, and we have as- 
sumed Q,m + f^A ~ 1 (used in the Millennium simulation) 
and matter-dominated era in computing dz/dt. As shown 
in Sec. 3, the parameters P and 7 in equation ([7| gener- 
ally depend on the halo mass. We find, however, that the 
mass dependence follows a simple power law independent of 
the redshift, and the simple analytic form in equation ([7| 
provides an excellent approximation for the mean mass ac- 
cretion rate as a function of redshift and halo mass: 



M) = 42Moyr~\ 



M 



(8) 



For completeness, the best fit for the medtan growth rate 
computed in the Millennium simulation is 



(*>. 



= 24.1A/Qyr"^ 



M 



1012 A//0 



We note that the overall amplitude of the mean is higher 
than the median due to the long, positive, M tail (see Fig. 
§. 

Fig.[5]compares the mean accretion rates of dark matter 
in Mq per year computed from the Millennium simulation 
(solid curves) and this formula (dashed curves) for haloes of 
mass 10^^A//o to lO^^A^o over the redshift range of and 5. 
The overall trend of the accretion rate is such that M /M has 
a weak dependence on M (oc Af"'^^^), and its dependence 
on redshift is approximately (1 -I- z)^'^" at low z and (1 + 
z)'^'^ at 2; > 1. This z-dependence is motivated by our 2- 
parameter form for M{z) and is more accurate than the 



c(l + 1.752)V^^™(l + z)^-hnA. 



(9) 



simple power law used in Genel et al. (20081, Neistein et al 
( 12006[ ), and [Neistein fc Dekel| ( |2008| ; our z ~ value, on 
the other hand, is consistent with theirs to within 20%. We 
have also computed M from the fitting form for the median 
MAH in the recent preprint by Z08. We found their Af to 
have a slightly steeper z-dependence than our equation [9] 
where their median value is within 20% of our median M at 
~ but exceeds ours by a factor of ~ 2 at 2 ~ 4. 

Along the right side of the vertical axis of Fig. [5] we la- 
bel the corresponding mean accretion rates of baryons, Mb, 
assuming a cosmic baryon-to-dark matter ratio of Q.b/Q.m ~ 
1/6. The results shown should be a reasonable approxima- 
tion for the mean rate of baryon mass that is being ac- 
creted at the virial radius of a dark matter halo of a given 
mass. Fig. |5] and equation (|8| indicate that this rate is 
Mh ~ 7Af0yr~^ for 10^^ Mq haloes today, and it increases 
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Figure 5. Mean mass accretion rate of dark matter as a function 
of rcdshift at halo mass lO", 1.5 X lO^^^ ^nd lO^^'Mo. The 
solid curves are computed from the Millennium haloes (except 
IQl^M©, which falls below our resolution limit of 1000 particles 
per halo) at a given mass (±20% range); the dashed curves show 
the accurate fit provided by eq. (jsj . The right side of the vertical 
axis labels the mean accretion rate of baryons, Mf,i assuming a 
cosmic baryon-to-dark matter ratio of ~ 1/6. The slight dip in M 
at 2 = is due to the artificial edge effect inherent in the stitch-3 
algorithm used to process the FOF merger trees. 



to 27, 69, and 140 A/eyr"^ for lO^^M© haloes at z = 1, 2, 
and 3, respectively. Since the infalling baryons are a reser- 
voir for the gas that fuels star formation, it is interesting to 
compare Mh with the mean star formation rates of different 
types of galaxies, e.g., M* ~ 4M0yr~^ for the Milky Way 
(e.g., [Diehl et al.||2006| , suggesting that about half of the 
infalling ~ 7 Mqyt~^ for Galactic-size haloes needs to 
be converted into stars. The relations among these different 
accretion rates and the implications will be investigated in 
a subsequent work. 

Having determined the mean rates, we show their dis- 
tributions and dispersions in Fig. [6] Four redshifts, z = 0, 
0.5, 1, and 2, and four ranges of halo masses (left to right 
panel) are shown. Both the differential (top panels) and cu- 
mulative (bottom panels) distributions of Mj, are plotted for 
comparison. Within each panel, the distribution of Mb at a 
given halo mass is seen to broaden significantly with increas- 
ing redshift. For instance, the (comoving) number density 
of 1.5 X IO^^Mq haloes with Mh > 250Mq yr"^ increases 
dramatically from 5 x 10"'' Mpc"^ at z = to 5 x 10"^ 
Mpc~^ at z = 2. At a given redshift, the distribution of 
Mb also broadens with increasing halo mass, although the 
distribution (and dispersion) of the ratio Mb /Mb is largely 
independent of mass. The latter is similar to the weak mass 
dependence of the mean M /M given by equation (Isl. 



5 FORMATION REDSHIFTS: MEAN AND 
DISPERSION 

It is well established that on average, more massive haloes 
form later than less massive haloes in the ACDM cosmol- 
ogy. The Millennium database provides sufficient statistics 
for us to quantify the distributions of the formation red- 
shift Zf and its mean and scatter over a wide range of halo 
masses (~ 10^^ to ~ lO^^M©). The formation redshift, along 



i 



10 .'^ 



with the late-time growth rate /3 — 7, can be thought of as 
two physically motivated quantities parameterizing the halo 
MAH. 

The distributions of for each type of MAHs for three 
halo mass bins are plotted in Fig. [7] The overall trend of de- 
creasing Zf with increasing halo mass is evident across the 
three panels. Within each panel, a correlation between Z[ 
and the MAH type is clearly seen. Nearly all haloes in the 
smallest few 2/ bins are Type II. This means that despite 
the fact that Type II dominates the highest mass bins, the 
haloes that constitute Type II are not merely the especially 
massive haloes which formed late, but also include less mas- 
sive haloes which formed late. On average, a Type II halo 
has a formation redshift 0.5 smaller than a typical halo. To 
a lesser degree. Type HI and Type IV are also distinct from 
the overall distribution. Both tend to form early. Type III 
more so than Type IV, and together the two types account 
for most of the haloes that formed early. 

Fig. |8] shows the mean formation redshift 2/ as a func- 
tion of halo mass for all 478,781 2 = Millennium haloes 
(leftmost panel) as well as for each type of MAH. As the 
scatter about the line is, to a good approximation, Gaussian, 
the 1(7 range about the line is also provided in the plots in 
each panel (hght shaded areas). From these shaded areas, 
it is clear that there is considerable scatter for the overall 
distribution. The relationship between Mo and 2/ is differ- 
ent from the overall distribution for all types except Type 
I, which suggests that the types discriminate by formation 
redshift to some extent. Also note that the separation of 
haloes into types also produces more limited scatter about 
the mean. 

To approximate the mass dependence of the mean and 
scatter of 2/, we use the linear form 



{2/) = alogio 



Mo 



Mo 



10i2Mp 



+ h, a-^, = clog„, +d 

(10) 

and find it to fit the simulation data accurately. Table [S] 
lists the best-fit coefficients for all the halo MAHs (above 
1000 particles at 2 = 0) and for each of the four types of 
MAHs shown in Fig. |8] Table [3] also includes the same fit 
performed for the fit parameters (/3— 7). The mean formation 
redshifts differ significantly among the types, with (zf) ~ 
0.6, 1.3, and 1.5 for Type II, I, and HI (plus IV), respectively, 
for galaxy-size haloes. The dependence of (zf) on mass is 
noticeably weak for Type II; the other types show similar 
mass dependence, where d{zf) /dlogM ranges from —0.23 
to -0.25. 

With the relationships between formation redshift and 
mass for each type, we can look at how these dependences 
relate to the basic halo characteristics given in Table [l] Re- 
call that Type II haloes were marked by steep growth at late 
times, which is captured by the very negative value of /3 — 7. 
Type HI haloes, on the other hand, have small values for 
/3 — 7, and thus grow slowly at late times. The relationships 
shown in Fig. [S] are then no surprise. Type II haloes are also 
associated with late formation times, while Type HI haloes 
tend to have formed quite early. 
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Figure 7. Distributions of the formation redshift Zf for three mass bins (left to right). Within each mass bin, the zj distribution for all 
haloes is plotted (solid black), as well as the distribution for each halo type. While the relative amplitudes of the distributions do change 
from one mass bin to another, the overall shapes remain similar across all masses. 
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Figure 8. Mass dependence of the mean (solid curves) and one standard deviation scatter (shaded regions) of the formation redshifts 
of the Millennium haloes. More massive haloes on average form more recently, but the scatter is large. One exception is Type II haloes 
that have a mean Zf of 0.5 independent of mass. Fits to the mean and scatter of Zf as a function of mass arc given in Tablejs] 



CORRELATIONS WITH HALO 
ENVIRONMENTS AND MAJOR MERGER 
FREQUENCIES 



Thus far we have discussed how the halo MAHs and mass 
accretion rates vary with halo mass and redshift. We have 
also shown that the mean 2/ depends on halo mass strongly, 
but the scatter in Zf does not depend strongly on the MAH 
type nor halo mass. In this section, we investigate if the 
mean and scatter in z/ are correlated with quantities other 
than halo mass. In particular, we ask if the shapes of MAHs 
(1) differ systematically between underdense vs overdense 
regions, and (2) are correlated with the time and frequency 
of major mergers and mass brought in by these events during 
a halo's lifetime. 



6.1 Environment 

An extensive discussion and tests of halo environments can 
be found in Fakhouri & Ma (20091. Four definitions of a 



halo's local environment based on the local mass density 
centered at the halo were compared. Three of them were 
computed using the dark matter particles in a sphere of ra- 
dius R centered at a halo, either with or without the central 
region carved out; the fourth definition was computed using 
the masses of only the haloes rather than all the dark mat- 
ter. Here we use Sji-pop, computed by subtracting out the 
FOF mass M of the central halo within a sphere of radius 
R: 



:5r 



M 

Vrp„ 



(11) 



where Vr is the volume of a sphere of radius R, and 5r is 
the mean mass overdensity within R. This measure makes 
no assumption about the central halo's shape. By taking out 
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Zf vs. Mo 1 — vs. Mq 

(fx) (7-/3) "-(7-/3) 

Overall -0.24a; + 1.26 -O.llx + 0.58 -0.25a: - 0.29 0.14a; + 0.64 
Type I -0.25X + 1.30 -0.12a; + 0.51 -0.22a: - 0.45 0.03a; + 0.21 
Type II -0.05X + 0.62 -0.04a; + 0.34 -0.12a: -1.15 0.14a; + 0.60 
Type III -0.23X + 1.56 -0.08a; + 0.43 -0.03a: - 0.20 O.Ola; + 0.11 
Type IV -0.23a; + 1.56 -0.08a' + 0.43 0.03a' + 0.15 0.36 

Table 3. Linear fits for the mass dependence of the mean formation redshift Zf, mean 7-/3, and their respective lo" scatter about the 
mean, where x = logiQ(Af()/10"'"^A/Q). 




l + S 



Figure 9. Normalized distributions in 5 are plotted for all haloes, 
as well as each type of halo. The most prominent feature is the 
clear separation between Type IV and all other types. The center 
of the Type IV is well to the right of the other types, meaning that 
Type IV haloes are predominantly found in denser environments 
than any other type. 



the mass of the central halo itself, this density was shown 
to be a more robust measure of the environment outside of 
a halo's virial radius. Otherwise, the halo mass itself dom- 
inates the density centered at massive haloes (see Fig. 1 of 



Fakhouri & Ma 2009 1 , and it becomes difficult to distin- 



guish whether any correlations are due to the mass or the 
larger-scale environment in which the halo resides. 

Fig. |9] shows the distribution of the environmental den- 
sities (evaluated at z = 0) for haloes within each MAH 
type and the total population. The distribution for Type 
IV haloes is quite distinct from all other distributions and is 
offset towards higher densities. Since Type IV haloes experi- 
ence very little mass growth or even mass loss at late times, 
denser environments appear to impede mass accretion onto 
haloes. 



6.2 Mass Growth due to Major Mergers 

Major mergers are more rare than minor mergers, but they 
can contribute to a significant fraction of a halo's final mass, 
and have a strong impact on halo structures and galaxy 
properties such as the star formation rate. 

A useful quantity for assessing the role of major mergers 
on halo MAHs is zimm, the redshift of the last major merger 
in a halo's history. Fig. |10| shows the fraction of haloes whose 
last major merger occurred at or before z. Type II haloes 
are seen to experience a major merger in the much more 
recent past than the other types: about 65% of them had 




Figure 10. Distribution of the redshift of the last merger major 
for the four types of haloes (colored curves) and slW z = haloes 
(solid black). The majority of the z = haloes with Type II 
MAH (dashed green curve) have experienced a major merger (of 
mass ratio > 0.33) very recently, whereas the last major merger 
occurred earlier than z = 1 for more than 75% of Type III and 
IV haloes. 



encountered a major merger within redshifts and 0.3, and 
only 25% of them had their last major merger before z = 1. 
In sharp contrast, only about 5% of Type III and IV haloes 
had a major merger later than z — 0.3, and over 75% of 
them had their last major merger earlier than z = 1. 

Another useful parameter for quantifying the role of ma- 
jor mergers in its MAH is F{£^ > ^min), which is the fraction 
of mass at zq that came from mergers above some progen- 
itor mass ratio ^min. We choose to define the mass ratio in 
relation to the mass of the progenitor at the time of merger, 
as opposed to being defined in relation to the halo's present 
mass. The exact value of F is strongly dependent upon the 
choice of ^min. The overall features, however, do not change 
significantly, so different values for ^min only change the val- 
ues for f (5 > Cmin), but leave the overall characteristics in 
place. 

Fig. [iT] shows the differential (top) and cumulative 
(bottom) distribution of the major merger mass fraction, 

> 0.33) for each type of MAHs. Type II haloes (dashed 
curves), which feature steep growth at late times, are seen 
to have the highest F among the four types. The distribu- 
tion peaks at F ^ 0.4, indicating that ~ 40% of their final 
mass was acquired through major mergers. Type I haloes 
(solid cyan curves), by contrast, feature a dearth of ma- 
jor mergers, which is unsurprising given the fact that large 
mergers are poorly handled by the simple exponential. Like- 
wise, Type IV haloes (dotted curves), which grow quickly 
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> 0.33) 

Figure 11. Differential (top) and cumulative (bottom) distribu- 
tions of the fraction of halo mass gained from major mergers of 
mass ratio ^ > 0.33. Much like the zj distributions, each Type is 
quite distinct from the overall distribution. Type I and III haloes 
feature few major mergers, while Type II and Type IV haloes tend 
to be dominated by major mergers. Type IV mergers also feature 
a noticeable rise in the highest bin, due to having gained more 
mass via major mergers across their history than they currently 
have. 



early on only to lose mass at late times, are dominated by 
major mergers with a number of Type IV haloes even hav- 
ing F{£^ > ^min) > 1. This occurs when a halo has less mass 
presently than it has gained overall via major merger. Type 
III haloes (dotted dashed curves) tend to be relatively ma- 
jor merger free, which is again to be expected due to the 
decelerating growth of Type III haloes at later times during 
which much of the overall mass is accreted. 



haloes increases sharply with increasing z at a given halo 
mass (see Fig. |6|. The mean halo formation redshift as a 
function of mass is given by equation ( 10 1 and Fig. [s] 

To facilitate the analysis of the halo MAH, we have clas- 
sified M{z) into four types based on their shapes. We have 
shown that only 20 to 30% of the Millennium haloes follow 
an exponential form ( "Type I" ) in their mass accretion his- 
tory M{z). Only one parameter is needed to specify their 
MAH, e.g., the formation redshift Zf. The formation red- 
shift depends strongly on halo mass, as expected for hierar- 
chical cosmological models such as the ACDM. The median 
Zf ranges from 1.3 for IO^^Mq haloes to 0.6 for 
haloes, and is dispersed over a range roughly equal to the 
median value for all masses. 

About 20% of galaxy-size and 60% of cluster-size haloes 
have late-time growth that is steeper than an exponential 
("Type H"). These haloes are formed more recently, with 
a median 2/ of about 0.5 for all masses. The redshift at 
which they experience the last major merger is also signifi- 
cantly later than the exponential haloes: about 50% of them 
have had the last major merger between z = and 0.3, as 
opposed to 10% of the rest of the haloes, including expo- 
nential haloes. Correspondingly, a higher fraction of Type 
II haloes' final mass is acquired through major mergers, e.g. 
60% of these haloes obtained more than 30% of their final 
mass from major mergers, whereas a little over 30% of all 
haloes obtained more than 30% of their final mass from ma- 
jor mergers, and fewer than 20% of exponential haloes did. 

The rest of the haloes have stunted late-time growth rel- 
ative to an exponential form. The median Zf ranges from 1.5 
at low mass to 0.8 at high mass. These haloes can be further 
separated into two groups (Type III and IV), where the two 
are primarily distinguished by the roles that major mergers 
have played in their growth; that is. Type III haloes tend to 
experience few major mergers, whereas Type IV haloes grew 
predominantly from major mergers at early redshifts. The 
MAHs of the two can be distinguished by the sharpness in 
the downturn of late time growth. Type IV haloes also live 
in somewhat denser environments, where the stronger tidal 
fields and more frequent interactions may have contributed 
to rapid accretion at early times followed by a slow down of 
their late time mass growth. 

Despite this diverse behavior of halo MAHs, we have 
found the individual M{z) to be well fit when a second pa- 
rameter is introduced (eq. |4|. To quantify the statistics of 
M{z), we have provided fits to the joint probability distri- 
bution of the two MAH parameters (3 and 7 in the Ap- 
pendix. These can be used to generate realizations of halo 
mass growth histories in semi-analytic models of galaxy for- 
mation that incorporate realistic scatters about the mean 
trends. 



7 CONCLUSIONS 

We have examined the mass growth histories of ~ 500, 000 
z = haloes and their progenitors in the Millennium simu- 
lation. The two-parameter function in equation Q provides 
a reasonable fit for the MAHs of these haloes, as shown in 
Fig. [2] The mean mass accretion rate of dark matter (and 
baryons) as a function of halo mass and redshift is well ap- 
proximated by equation (|8|, as shown in Fig. [5] The distri- 
butions of M are broad, and the number density of high-M 
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APPENDIX A: JOINT DISTRIBUTION OF /3 
AND 7 

We have seen that halo MAHs are well-fit by equation Q 
with two parameters /? and 7. Applying this fit to haloes 
in the Millennium simulation yields a joint distribution ol 
j3 and 7. In this appendix we provide a fitting form to this 
distribution as a function of /3, 7, and halo mass that is 
intended to allow the reader to generate rapidly a mock cat- 
alog of MAH tracks. We find that a straightforward rejec- 
tion method can generate 300,000 mock MAH trajectories 
in under a minute on a standard laptop. The mean proper- 
ties of the resulting trajectories match the mean properties 
of the Millennium trajectories at the 10% level. The fitting 
forms presented below are chosen for the practical purpose 
of matching the underlying distribution as closely as possi- 
ble. 

We find that 95.34% of the haloes occupy a smooth re- 
gion in the {(3, 7) plane shown in the left panel of Fig. |A1[ 
The remaining 4.66% of the haloes live along a distinct line 



with 7 = and /3 ^ 0, where the distribution of (3 is shown 
in the right panel of Fig. |A1| That is, their MAHs are bet- 
ter approximated by a power law in 1 4- 2 rather than an 
exponential. Interestingly, this 95.34% vs 4.66% division is 
independent of halo mass, even though the shape of the 
distributions generally depends on mass. For accuracy, we 
choose to separate the distribution of /3 and 7 into these 
two components and fit to them separately. 

For the 4.66% of haloes with 7 = 0, their /3 distribution 
is well approximated by 

dP _x- 

where 

0. 6335/3+0. 2626Af°- 1^3^ 



(Al) 



X = 7.443 c" 



2.852M-°'°^''^^ (A2) 
^^Mq . This fit is valid in the range - 10 < 



and M = Mhaio/10 
0. 

For the 95.34% of haloes with 7 > 10 



, the joint dis- 



tribution in (3 and 7 is well approximated by 



where 
X = 

Y = 
Y' = 

Ti = 
T2 = 



dfddj 

(-1.722 - 0.1568/3 + 0.007592/3^)(l - T2) + 

(1.242 + 0.3138/3 - 0.01399/3^)T2 

13.39(1-1. 224 tanh(1.043F')](l-0. 08018/3) 

7 - (28.85 + 0.4537/3)(l - Ti) - 

(28.38 + 0.7624/3)ri + 29.21M~°°°^^^^ 

0.5[1 + tanh(l. 174/3)] 

0.5[1 + tanh(0.7671/3 - 0.1269)] . 



(A3) 



This is valid in the range -8 < ^ < 12, < V" < 3. Since 
the rejection method does not require a normalized PDF for 
input, we leave these probability distributions unnormalized. 

Fig. |A2| illustrates that Monte Carlo realizations gen- 
erated from the probability distributions above (dotted 
curves) reproduce accurately the formation redshift distri- 
butions (left panel) and the mean MAHs (right panel) ob- 
tained from the (/3, 7) fits to the Millennium MAHs (dashed 
curves), and both match closely the results computed di- 
rectly from the Millennium simulation (solid curves). 
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Figure Al. Probability distributions of (3 and 7 for haloes of mass 1.5 X lO^^M©. The left panel is the PDF of /3 and 7 for haloes with 
7 > 10-3 (95.34% of all haloes). The shaded 2D histo gram presents the distribution of fi and 7 obtained from Millennium (see color scale 
for units) with corresponding contours drawn at the 0.005, 0.01,0.1 and 1 levels (black to white). The appropriately normalized fitting 
form D2D is overlaid as red dashed contour lines. The right panel is the PDF of (3 for haloes with 7 = (4.66%). The appropriately 
normalized fitting form Did is overlaid in red. For both panels, the background shading corresponds to regions of phase space denote 
type I (grey), type II (red), type III (green) and type IV (blue) MAHs. 




Figure A2. Distribution of the formation redshifts (left) and the mean MAHs (right) from the simulation directly (solid), the (/3,7) 
fits to the simulation MAHs (dashed), and a Monte Carlo ensemble of 300,000 halos per mass bin generated using equations (A1)-(A3) 
(dotted). In each panel, three mass bins are shownrlO^'^M© (blue), 10^"^ Mq (green), and IO^^Mq (red). Note that the dashed and dotted 
curves are almost indistinguishable in the right panel. 
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